Information transfer and nontrivial collective behavior 
in chaotic coupled map networks 
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The emergence of nontrivial collective behavior in networks of coupled chaotic maps is investigated 
by means of a nonlinear mutual prediction method. The resulting prediction error is used to measure 
the amount of information that a local unit possesses about the collective dynamics. Applications 
to locally and globally coupled map systems are considered. The prediction error exhibits phase 
transitions at critical values of the coupling for the onset of ordered collective behavior in these 
networks. This information measure may be used as an order parameter for the characterization of 
complex behavior in extended chaotic systems. 

PACS numbers: 05.45.-a, 02.50.-r 
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Much interest has recently been directed to under- 
standing the phenomenon of emergence of nontrivial col- 
lective behavior in systems of interacting chaotic ele- 
ments Nontrivial collective behavior (NTCB) is 
characterized by a well-defined evolution of macroscopic 
quantities emerging out of local chaos. Models based 
on coupled map networks have been widely used in the 
investigation of cooperative phenomena that appear in 
many extended chaotic dynamical systems g{ . In partic- 
ular, NTCB has been studied in coupled maps on regular 
Euclidean lattices (lj g|, in one-dimensional lattices j|, 
fractal geometries ||, and globally coupled map systems 
H 0> H> a- ^ n tms ar ticl e > we investigate the information 
transfer between the local and global levels of coupled 
map network models as a condition for self-organization 
in spatiotemporal systems. One may ask the question: 
how much information does a local unit possess about 
the collective dynamics of a system?; or, how does the 
information flow depend on parameters of the system?. 

To analyze how global order in spatiotemporal sys- 
tems can arise out of local chaos, we consider a system 
of N interacting elements, where the state of element 
i (i = 1,2, ... ,N) at discrete time t is denoted by x\. 
The evolution of each element is assumed to depend on 
its own local map dynamics and on its interaction with 
other elements in the network, where the strength of the 
interactions is given by a coupling parameter. The collec- 
tive dynamics of the system at time t may be described 
by some statistical variable h t , such as the mean field. 
The information transfer between local (x\) and macro- 
scopic (h t ) variables is analyzed by comparing their cor- 
responding time evolutions. We adopt a simple compu- 
tational technique based on a mutual nonlinear modeling 
Jl0| , O. This method makes use of the temporal evolu- 
tion of a driven variable (which is receiving information) 



to infer characteristics of the driver variable. In our case, 
for different values of the coupling parameter we record 
time series containing the simultaneous evolution of the 
macroscopic quantity h t and of a chosen local variable 
x\. 

The trajectory of the series x\ is reconstructed in 
an embedding space of dimension d as a collection of 



vectors (x\,x 
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Then, for each of these 



vectors, we systematically look for its nearest neigh- 
bor (Xp, x'' p _ 1 , . . . ,x l p _ d+1 ) in the Euclidean distance, as 
shown in Fig. 1. The root-mean-square error over h t at 
the embedding dimension d is computed as 
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where h p is the value of the macroscopic variable that 
bears the time index of xt, and a is the standard devia- 
tion of ht- 

The prediction error given by Eq. (|l|) can be inter- 
preted as a measure of the information that the series 
x\ possesses about the macroscopic variable h t . In this 
context, small values of Ed(h\x) imply that the dynamics 
of the global quantity h t is embedded in the local evolu- 
tion. Larger mutual prediction errors indicate that the 
two time series become more independent of each other. 
In particular, for two totally independent random series 
the error Eq. (Q) has a value of -\/2. 

As applications, we have calculated the error E<i(h\x) 
as a function of the coupling strength in coupled map 
systems exhibiting NTCB. This quantity was computed 
for several embedding dimensions d and the curve giving 
minimal errors was selected in each case. 

The first example is a system of maps subjected to 
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FIG. 1: (a) Bifurcation diagram h t vs. e for homogeneous 
globally coupled maps, with b = 0. Four different phases are 
observed as e is varied: turbulent (T), periodic (P), chaotic 
bands (C), and synchronized (S). (b) Ez(h\x) vs e for this sys- 
tem, (c) Bifurcation diagram of a local mapa;J vs. e, exposing 
the underlying chaotic dynamics. 
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FIG. 2: (a) Bifurcation diagram h t vs. e for homogeneous 
globally coupled maps, with b = 0. Four different phases are 
observed as e is varied: turbulent (T), periodic (P), chaotic 
bands (C), and synchronized (S). (b) Es(h\x) vs e for this 
system, (c) Bifurcation diagram of a local map x\ vs. e, 
exposing the underlying chaotic dynamics. 



global coupling defined as 

N 

i=i 

where the function fi(x\) describes the local dynamics of 
element i, and e is the coupling parameter. The usual ho- 
mogeneous globally coupled map system fl corresponds 
to having the same local function for all the elements, 
i.e., fi{x\) — f(x\). As local dynamics, we employ the 
logarithmic map f(x) = b + In |a;| [Q, where 6 is a real 
parameter. This map does not belong to the standard 
classes of universality of unimodal or bounded maps. Ro- 
bust chaos occurs in the parameter interval b G [—1,1], 
with no periodic windows and no separated chaotic bands 
on this interval fl2]| . In a second example, heterogeneity 
in the local maps in Eq. (Jfy is introduced by taking 
fi{x\) = bi + In \ x\\, with the values of the parameters bi 
distributed in [—1, 1]. 

The last example is a one-dimensional, homogeneous, 
diffusively coupled logarithmic map lattice given by 

x\ +1 = (1 - e)f{x\) + \ [f{x\ +1 ) + f{x\- 1 )] , (3) 

where periodic boundary conditions are assumed in Eq. 
(§■ 

As the macroscopic variable for these systems we con- 
sider the instantaneous mean field, defined as 

1 N 



Each of the above examples presents nontrivial collec- 
tive behavior in some range of their parameters. Figure 
2(a) shows the bifurcation diagram of the mean field h t of 
the homogeneous globally coupled map system, Eq. (||), 
as a function of the coupling strength e ||] . The local pa- 
rameter is fixed at b — for all maps and the system size 
is N — 10 4 . For each value of e, the mean field was cal- 
culated at each time step during a run starting from ran- 
dom initial conditions on the local maps, uniformly dis- 
tributed on the interval [—8,4], after discarding the tran- 
sients. When the local parameter b is in the range [—1,1], 
the elements x\ are chaotic and desynchronized (see Fig. 
2(c)). However, the mean field in Fig. 2(a) reveals the 
existence of global periodic attractors for some intervals 
of the coupling. Different collective states emerge as a 
function of the coupling e: a turbulent phase, where ht 
manifests itself as a fixed point, following the standard 
statistical behavior of uncorrelated disordered variables; 
collective periodic states; collective chaotic bands; and 
chaotic synchronization In this representation, col- 
lective periodic states at a given value of the coupling e 
appear as sets of vertical segments which correspond to 
intrinsic fluctuations of the periodic orbits of the mean 
field. Increasing the system size N does not decrease 
the amplitude of the collective periodic orbits. More- 
over, when N is increased the width of the segments that 
make a periodic orbit in the bifurcation diagrams such 
as in Fig. 2(a) shrink, indicating that the global peri- 
odic attractors become better defined in the large sys- 
tem limit. Figure 2(b) shows E^{h\x) vs. e. An abrupt 
change in the value of the quantity E^(h\x) can clearly 
be seen at a critical value of the coupling e c ~ 0.21, a be- 
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havior characteristic of a first order phase transition. At 
this critical value of the coupling the collective behavior 
changes from a turbulent phase to a periodic collective 
state, as observed in Fig. 2(a). The error for e < e c is 
large, indicating that there is no appreciable information 
sharing between the local and macroscopic levels when 
the system is in the turbulent phase. For e > e c the er- 
ror drops discontinuously to very small values. The local 
unit suddenly becomes "aware" of the collective dynam- 
ics; the time series of a single map is good enough to 
provide assertive predictions of the mean field evolution. 
Thus, there is a large amount of information transfer 
from the collective dynamics to each of the elements in 
the network, even before synchronization is achieved. In- 
creasing the coupling beyond the synchronization region 
leads again to a turbulent state of the system and to a 
large value of the error E 3 (h\x). 
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FIG. 3: (a) Bifurcation diagram ht vs. e for heterogeneous 
globally coupled maps; system size is N = 10 4 . (b) E z {h\x) 
vs e for this network. 

Figure 3(a) shows the bifurcation diagram of h t vs. 
e for the globally coupled heterogeneous map lattice. 
In this case the local parameters hi are set at random 
with a uniform distribution in the chaotic interval, i.e., 
hi e [—1,1]. Again the local dynamics are chaotic, yet 
collective periodic behavior arise in some windows of the 
coupling parameter. Figure 3(b) shows the error E^{h\x) 
vs. e for the heterogeneous globally coupled system. In 
this case the prediction error stays large up to a critical 
value of the coupling e c ~ 0.04 and then decreases con- 
tinuously for e > e c . The decrease in the prediction error 
resembles a second order phase transition. The emer- 
gence of collective periodic behavior is manifested in the 
low values of the error for e > e c . The disorder introduced 
in the globally coupled network by the local heterogene- 
ity can be detected by the quantity E^{h\x) as a change 
in the character of the transition to collective behavior 
when compared to a similar transition in a homogeneous 
globally coupled system, Fig. 2(a). 

The emergence of nontrivial collective behavior can 



also be observed in the one-dimensional coupled map lat- 
tice, Eq. (||), as shown in Fig. 4(a). In this case the 
system presents only a turbulent (statistical fixed point) 
phase and a period-two collective state 0. Figure 4(b) 
shows E2{h\x) as a function of the coupling in the one- 
dimensional coupled map lattice. There is again a de- 
crease in the error at the critical value of the coupling 
strength for which the transition from turbulence to pe- 
riodic collective states occurs. We have observed simi- 
lar transitions from large to small values in the quantity 
Ed(h\x) at the onset of NTCB in networks having other 
connecting topologies, as well as when employing uni- 
modal chaotic maps as local dynamics. 
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FIG. 4: (a) Bifurcation diagram ht vs. e for one- dimensional 
coupled map lattice; N = 10 4 and 6 = 0. (b) E2(h\x) vs e for 
this lattice. 

The logarithmic map has been employed as local dy- 
namics in the above examples because the emergence of 
ordered collective behavior in those coupled systems can- 
not be attributed to the existence of windows of period- 
icity nor to chaotic band splitting in the local dynamics. 
These systems can be chaotic at a local level and simul- 
taneously periodic at a macroscopic level. Thus, there 
should exist some global information sharing among the 
elements of the networks that leads to a collective organi- 
zation besides the trivial synchronization. The informa- 
tion transfer required for the onset of nontrivial collective 
behavior takes place at some specific values of the pa- 
rameters of the systems. The observed decreasing of the 
errors at the transition to nontrivial collective behavior 
can be interpreted as a manifestation of the emergence 
of organization in these systems. It should be noticed 
that in all cases the dynamics of the elements do not ex- 
perience notable change before and after the transition 
to nontrivial collective behavior, since local dynamics is 
always chaotic. 

In conclusion, we have shown that the error Ed(h\x) 
is a useful quantity to characterize the transition to or- 
dered collective behavior in chaotic spatiotemporal sys- 
tems. Connectivity and coupling strengths are the mech- 
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anism for information flow in networks of dynamical 
units. However, our results suggest that transference of 
the information that is relevant for the emergence of col- 
lective organization in systems of interacting chaotic ele- 
ments is associated to low values of Ed{h\x). Finally, the 
exploration of a possible relationship between Ed{h\x) 
and other quantities used to study coupled chaotic map 
systems, such as the collective Lyapunov exponent || or 
transfer entropy p| , is an interesting problem for future 
research. 
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